An IL-1ß-driven neutrophil-stromal cell axis fosters a BAFF-rich protumor microenvironment in individuals with multiple myeloma.

Human bone marrow permanently harbors high numbers of neutrophils, and a tumor-supportive bias of these cells could significantly impact bone marrow-confined malignancies. In individuals with multiple myeloma, the bone marrow is characterized by inflammatory stromal cells with the potential to influence neutrophils. We investigated myeloma-associated alterations in human marrow neutrophils and the impact of stromal inflammation on neutrophil function. Mature neutrophils in myeloma marrow are activated and tumor supportive and transcribe increased levels of IL1B and myeloma cell survival factor TNFSF13B (BAFF). Interactions with inflammatory stromal cells induce neutrophil activation, including BAFF secretion, in a STAT3-dependent manner, and once activated, neutrophils gain the ability to reciprocally induce stromal activation. After first-line myeloid-depleting antimyeloma treatment, human bone marrow retains residual stromal inflammation, and newly formed neutrophils are reactivated. Combined, we identify a neutrophil-stromal cell feed-forward loop driving tumor-supportive inflammation that persists after treatment and warrants novel strategies to target both stromal and immune microenvironments in multiple myeloma.

A Single-Cell Taxonomy Predicts Inflammatory Niche Remodeling to Drive Tissue Failure and Outcome in Human AML.

Cancer initiation is orchestrated by an interplay between tumor-initiating cells and their stromal/immune environment. Here, by adapted single-cell RNA sequencing, we decipher the predicted signaling between tissue-resident hematopoietic stem/progenitor cells (HSPC) and their neoplastic counterparts with their native niches in the human bone marrow. LEPR+ stromal cells are identified as central regulators of hematopoiesis through predicted interactions with all cells in the marrow. Inflammatory niche remodeling and the resulting deprivation of critical HSPC regulatory factors are predicted to repress high-output hematopoietic stem cell subsets in NPM1-mutated acute myeloid leukemia (AML), with relative resistance of clonal cells. Stromal gene signatures reflective of niche remodeling are associated with reduced relapse rates and favorable outcomes after chemotherapy across all genetic risk categories. Elucidation of the intercellular signaling defining human AML, thus, predicts that inflammatory remodeling of stem cell niches drives tissue repression and clonal selection but may pose a vulnerability for relapse-initiating cells in the context of chemotherapeutic treatment. SIGNIFICANCE: Tumor-promoting inflammation is considered an enabling characteristic of tumorigenesis, but mechanisms remain incompletely understood. By deciphering the predicted signaling between tissue-resident stem cells and their neoplastic counterparts with their environment, we identify inflammatory remodeling of stromal niches as a determinant of normal tissue repression and clinical outcomes in human AML. See related commentary by Lisi-Vega and Méndez-Ferrer, p. 349. This article is featured in Selected Articles from This Issue, p. 337.

Hematopoietic Cell Autonomous Disruption of Hematopoiesis in a Germline Loss-of-function Mouse Model of RUNX1-FPD.

RUNX1 familial platelet disorder (RUNX1-FPD) is a hematopoietic disorder caused by germline loss-of-function mutations in the RUNX1 gene and characterized by thrombocytopathy, thrombocytopenia, and an increased risk of developing hematologic malignancies, mostly of myeloid origin. Disease pathophysiology has remained incompletely understood, in part because of a shortage of in vivo models recapitulating the germline RUNX1 loss of function found in humans, precluding the study of potential contributions of non-hematopoietic cells to disease pathogenesis. Here, we studied mice harboring a germline hypomorphic mutation of one Runx1 allele with a loss-of-function mutation in the other Runx1 allele (Runx1 L148A/- mice), which display many hematologic characteristics found in human RUNX1-FPD patients. Runx1 L148A/- mice displayed robust and pronounced thrombocytopenia and myeloid-biased hematopoiesis, associated with an HSC intrinsic reconstitution defect in lymphopoiesis and expansion of myeloid progenitor cell pools. We demonstrate that specific deletion of Runx1 from bone marrow stromal cells in Prrx1-cre;Runx1 fl/fl mice did not recapitulate these abnormalities, indicating that the hematopoietic abnormalities are intrinsic to the hematopoietic lineage, and arguing against a driving role of the bone marrow microenvironment. In conclusion, we report a RUNX1-FPD mouse model faithfully recapitulating key characteristics of human disease. Findings do not support a driving role of ancillary, non-hematopoietic cells in the disruption of hematopoiesis under homeostatic conditions.

Gata2-regulated Gfi1b expression controls endothelial programming during endothelial-to-hematopoietic transition.

The first hematopoietic stem cells (HSCs) are formed through endothelial-to-hematopoietic transition (EHT) during embryonic development. The transcription factor GATA2 is a crucial regulator of EHT and HSC function throughout life. Because patients with GATA2 haploinsufficiency have inborn mutations, prenatal defects are likely to influence disease development. In mice, Gata2 haploinsufficiency (Gata2+/-) reduces the number and functionality of embryonic hematopoietic stem and progenitor cells (HSPCs) generated through EHT. However, the embryonic HSPC pool is heterogeneous and the mechanisms underlying this defect in Gata2+/- embryos remain unclear. Here, we investigated whether Gata2 haploinsufficiency selectively affects a cellular subset undergoing EHT. We showed that Gata2+/- HSPCs initiate, but cannot fully activate, hematopoietic programming during EHT. In addition, due to the reduced activity of the endothelial repressor Gfi1b, Gata2+/- HSPCs cannot repress endothelial identity to complete maturation. Finally, we showed that hematopoietic-specific induction of gfi1b could restore HSC production in gata2b-null (gata2b-/-) zebrafish embryos. This study illustrates the pivotal role of Gata2 in the regulation of the transcriptional network governing HSPC identity throughout the EHT.

APOBEC mutagenesis is a common process in normal human small intestine.

APOBEC mutational signatures SBS2 and SBS13 are common in many human cancer types. However, there is an incomplete understanding of its stimulus, when it occurs in the progression from normal to cancer cell and the APOBEC enzymes responsible. Here we whole-genome sequenced 342 microdissected normal epithelial crypts from the small intestines of 39 individuals and found that SBS2/SBS13 mutations were present in 17% of crypts, more frequent than most other normal tissues. Crypts with SBS2/SBS13 often had immediate crypt neighbors without SBS2/SBS13, suggesting that the underlying cause of SBS2/SBS13 is cell-intrinsic. APOBEC mutagenesis occurred in an episodic manner throughout the human lifespan, including in young children. APOBEC1 mRNA levels were very high in the small intestine epithelium, but low in the large intestine epithelium and other tissues. The results suggest that the high levels of SBS2/SBS13 in the small intestine are collateral damage from APOBEC1 fulfilling its physiological function of editing APOB mRNA.

Prognostic Value of FLT3-Internal Tandem Duplication Residual Disease in Acute Myeloid Leukemia.

PURPOSE: The applicability of FLT3-internal tandem duplications (FLT3-ITD) for assessing measurable residual disease (MRD) in acute myeloid leukemia (AML) in complete remission (CR) has been hampered by patient-specific duplications and potential instability of FLT3-ITD during relapse. Here, we comprehensively investigated the impact of next-generation sequencing (NGS)-based FLT3-ITD MRD detection on treatment outcome in a cohort of patients with newly diagnosed AML in relation to established prognostic factors at diagnosis and other MRD measurements, ie, mutant NPM1 and multiparameter flow cytometry. METHODS: In 161 patients with de novo FLT3-ITD AML, NGS was performed at diagnosis and in CR after intensive remission induction treatment. FLT3-ITD MRD status was correlated with the cumulative incidence of relapse and overall survival (OS). RESULTS: NGS-based FLT3-ITD MRD was present in 47 of 161 (29%) patients with AML. Presence of FLT3-ITD MRD was associated with increased risk of relapse (4-year cumulative incidence of relapse, 75% FLT3-ITD MRD v 33% no FLT3-ITD MRD; P < .001) and inferior OS (4-year OS, 31% FLT3-ITD MRD v 57% no FLT3-ITD MRD; P < .001). In multivariate analysis, detection of FLT3-ITD MRD in CR confers independent prognostic significance for relapse (hazard ratio, 3.55; P < .001) and OS (hazard ratio 2.51; P = .002). Strikingly, FLT3-ITD MRD exceeds the prognostic value of most generally accepted clinical and molecular prognostic factors, including the FLT3-ITD allelic ratio at diagnosis and MRD assessment by NGS-based mutant NPM1 detection or multiparameter flow cytometry. CONCLUSION: NGS-based detection of FLT3-ITD MRD in CR identifies patients with AML with profound risk of relapse and death that outcompetes the significance of most established prognostic factors at diagnosis and during therapy, and furnishes support for FLT3-ITD as a clinically relevant biomarker for dynamic disease risk assessment in AML.

Myeloid cells promote interferon signaling-associated deterioration of the hematopoietic system.

Innate and adaptive immune cells participate in the homeostatic regulation of hematopoietic stem cells (HSCs). Here, we interrogate the contribution of myeloid cells, the most abundant cell type in the mammalian bone marrow, in a clinically relevant mouse model of neutropenia. Long-term genetic depletion of neutrophils and eosinophils results in activation of multipotent progenitors but preservation of HSCs. Depletion of myeloid cells abrogates HSC expansion, loss of serial repopulation and lymphoid reconstitution capacity and remodeling of HSC niches, features previously associated with hematopoietic aging. This is associated with mitigation of interferon signaling in both HSCs and their niches via reduction of NK cell number and activation. These data implicate myeloid cells in the functional decline of hematopoiesis, associated with activation of interferon signaling via a putative neutrophil-NK cell axis. Innate immunity may thus come at the cost of system deterioration through enhanced chronic inflammatory signaling to stem cells and their niches.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus9 (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Clonal diversification and histogenesis of malignant germ cell tumours.

Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells. They are uniquely able to recapitulate embryonic and extraembryonic tissues, which carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we study the relationship of histogenesis and clonal diversification in GCTs by analysing the genomes and transcriptomes of 547 microdissected histological units. We find no correlation between genomic and histological heterogeneity. However, we identify unifying features including the retention of fetal developmental transcripts across tissues, expression changes on chromosome 12p, and a conserved somatic evolutionary sequence of whole genome duplication followed by clonal diversification. While this pattern is preserved across all GCTs, the developmental timing of the duplication varies between prepubertal and postpubertal cases. In addition, tumours of younger children exhibit distinct substitution signatures which may lend themselves as potential biomarkers for risk stratification. Our findings portray the extensive diversification of GCT tissues and genetic subclones as randomly distributed, while identifying overarching transcriptional and genomic features.

Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Inherited MUTYH mutations cause elevated somatic mutation rates and distinctive mutational signatures in normal human cells.

Cellular DNA damage caused by reactive oxygen species is repaired by the base excision repair (BER) pathway which includes the DNA glycosylase MUTYH. Inherited biallelic MUTYH mutations cause predisposition to colorectal adenomas and carcinoma. However, the mechanistic progression from germline MUTYH mutations to MUTYH-Associated Polyposis (MAP) is incompletely understood. Here, we sequence normal tissue DNAs from 10 individuals with MAP. Somatic base substitution mutation rates in intestinal epithelial cells were elevated 2 to 4-fold in all individuals, except for one showing a 31-fold increase, and were also increased in other tissues. The increased mutation burdens were of multiple mutational signatures characterised by C > A changes. Different mutation rates and signatures between individuals are likely due to different MUTYH mutations or additional inherited mutations in other BER pathway genes. The elevated base substitution rate in normal cells likely accounts for the predisposition to neoplasia in MAP. Despite ubiquitously elevated mutation rates, individuals with MAP do not display overt evidence of premature ageing. Thus, accumulation of somatic mutations may not be sufficient to cause the global organismal functional decline of ageing.

Clonal evolution after treatment pressure in multiple myeloma: heterogenous genomic aberrations and transcriptomic convergence.

We investigated genomic and transcriptomic changes in paired tumor samples of 29 in-house multiple myeloma (MM) patients and 28 patients from the MMRF CoMMpass study before and after treatment. A change in clonal composition was found in 46/57 (82%) of patients, and single-nucleotide variants (SNVs) increased from median 67 to 86. The highest increase in prevalence of genetic aberrations was found in RAS genes (60% to 72%), amp1q21 (18% to 35%), and TP53 (9% to 18%). The SBS-MM1 mutation signature was detected both in patients receiving high and low dose melphalan. A total of 2589 genes were differentially expressed between early and late samples (FDR < 0.05). Gene set enrichment analysis (GSEA) showed increased expression of E2F, MYC, and glycolysis pathways and a decreased expression in TNF-NFkB and TGFbeta pathways in late compared to early stage. Single sample GSEA (ssGSEA) scores of differentially expressed pathways revealed that these changes were most evident in end-stage disease. Increased expression of several potentially targetable genes was found at late disease stages, including cancer-testis antigens, XPO1 and ABC transporters. Our study demonstrates a transcriptomic convergence of pathways supporting increased proliferation and metabolism during disease progression in MM.

Mutational landscape of normal epithelial cells in Lynch Syndrome patients.

Lynch Syndrome (LS) is an autosomal dominant disease conferring a high risk of colorectal cancer due to germline heterozygous mutations in a DNA mismatch repair (MMR) gene. Although cancers in LS patients show elevated somatic mutation burdens, information on mutation rates in normal tissues and understanding of the trajectory from normal to cancer cell is limited. Here we whole genome sequence 152 crypts from normal and neoplastic epithelial tissues from 10 LS patients. In normal tissues the repertoire of mutational processes and mutation rates is similar to that found in wild type individuals. A morphologically normal colonic crypt with an increased mutation burden and MMR deficiency-associated mutational signatures is identified, which may represent a very early stage of LS pathogenesis. Phylogenetic trees of tumour crypts indicate that the most recent ancestor cell of each tumour is already MMR deficient and has experienced multiple cycles of clonal evolution. This study demonstrates the genomic stability of epithelial cells with heterozygous germline MMR gene mutations and highlights important differences in the pathogenesis of LS from other colorectal cancer predisposition syndromes.

Somatic mutation rates scale with lifespan across mammals.

The rates and patterns of somatic mutation in normal tissues are largely unknown outside of humans1-7. Comparative analyses can shed light on the diversity of mutagenesis across species, and on long-standing hypotheses about the evolution of somatic mutation rates and their role in cancer and ageing. Here we performed whole-genome sequencing of 208 intestinal crypts from 56 individuals to study the landscape of somatic mutation across 16 mammalian species. We found that somatic mutagenesis was dominated by seemingly endogenous mutational processes in all species, including 5-methylcytosine deamination and oxidative damage. With some differences, mutational signatures in other species resembled those described in humans8, although the relative contribution of each signature varied across species. Notably, the somatic mutation rate per year varied greatly across species and exhibited a strong inverse relationship with species lifespan, with no other life-history trait studied showing a comparable association. Despite widely different life histories among the species we examined-including variation of around 30-fold in lifespan and around 40,000-fold in body mass-the somatic mutation burden at the end of lifespan varied only by a factor of around 3. These data unveil common mutational processes across mammals, and suggest that somatic mutation rates are evolutionarily constrained and may be a contributing factor in ageing.

Molecular characterization of mutant TP53 acute myeloid leukemia and high-risk myelodysplastic syndrome.

Substantial heterogeneity within mutant TP53 acute myeloid leukemia (AML) and myelodysplastic syndrome with excess of blast (MDS-EB) precludes the exact assessment of prognostic impact for individual patients. We performed in-depth clinical and molecular analysis of mutant TP53 AML and MDS-EB to dissect the molecular characteristics in detail and determine its impact on survival. We performed next-generation sequencing on 2200 AML/MDS-EB specimens and assessed the TP53 mutant allelic status (mono- or bi-allelic), the number of TP53 mutations, mutant TP53 clone size, concurrent mutations, cytogenetics, and mutant TP53 molecular minimal residual disease and studied the associations of these characteristics with overall survival. TP53 mutations were detected in 230 (10.5%) patients with AML/MDS-EB with a median variant allele frequency of 47%. Bi-allelic mutant TP53 status was observed in 174 (76%) patients. Multiple TP53 mutations were found in 49 (21%) patients. Concurrent mutations were detected in 113 (49%) patients. No significant difference in any of the aforementioned molecular characteristics of mutant TP53 was detected between AML and MDS-EB. Patients with mutant TP53 have a poor outcome (2-year overall survival, 12.8%); however, no survival difference between AML and MDS-EB was observed. Importantly, none of the molecular characteristics were significantly associated with survival in mutant TP53 AML/MDS-EB. In most patients, TP53 mutations remained detectable in complete remission by deep sequencing (73%). Detection of residual mutant TP53 was not associated with survival. Mutant TP53 AML and MDS-EB do not differ with respect to molecular characteristics and survival. Therefore, mutant TP53 AML/MDS-EB should be considered a distinct molecular disease entity.

Systematic Profiling of DNMT3A Variants Reveals Protein Instability Mediated by the DCAF8 E3 Ubiquitin Ligase Adaptor.

Clonal hematopoiesis is a prevalent age-related condition associated with a greatly increased risk of hematologic disease; mutations in DNA methyltransferase 3A (DNMT3A) are the most common driver of this state. DNMT3A variants occur across the gene with some particularly associated with malignancy, but the functional relevance and mechanisms of pathogenesis of the majority of mutations are unknown. Here, we systematically investigated the methyltransferase activity and protein stability of 253 disease-associated DNMT3A mutations, and found that 74% were loss-of-function mutations. Half of these variants exhibited reduced protein stability and, as a class, correlated with greater clonal expansion and acute myeloid leukemia development. We investigated the mechanisms underlying the instability using a CRISPR screen and uncovered regulated destruction of DNMT3A mediated by the DCAF8 E3 ubiquitin ligase adaptor. We establish a new paradigm to classify novel variants that has prognostic and potential therapeutic significance for patients with hematologic disease. SIGNIFICANCE: DNMT3A has emerged as the most important epigenetic regulator and tumor suppressor in the hematopoietic system. Our study represents a systematic and high-throughput method to characterize the molecular impact of DNMT3A missense mutations and the discovery of a regulated destruction mechanism of DNMT3A offering new prognostic and future therapeutic avenues.See related commentary by Ma and Will, p. 23.This article is highlighted in the In This Issue feature, p. 1.

Increased somatic mutation burdens in normal human cells due to defective DNA polymerases.

Mutation accumulation in somatic cells contributes to cancer development and is proposed as a cause of aging. DNA polymerases Pol ε and Pol δ replicate DNA during cell division. However, in some cancers, defective proofreading due to acquired POLE/POLD1 exonuclease domain mutations causes markedly elevated somatic mutation burdens with distinctive mutational signatures. Germline POLE/POLD1 mutations cause familial cancer predisposition. Here, we sequenced normal tissue and tumor DNA from individuals with germline POLE/POLD1 mutations. Increased mutation burdens with characteristic mutational signatures were found in normal adult somatic cell types, during early embryogenesis and in sperm. Thus human physiology can tolerate ubiquitously elevated mutation burdens. Except for increased cancer risk, individuals with germline POLE/POLD1 mutations do not exhibit overt features of premature aging. These results do not support a model in which all features of aging are attributable to widespread cell malfunction directly resulting from somatic mutation burdens accrued during life.

Convergent somatic mutations in metabolism genes in chronic liver disease.

The progression of chronic liver disease to hepatocellular carcinoma is caused by the acquisition of somatic mutations that affect 20-30 cancer genes1-8. Burdens of somatic mutations are higher and clonal expansions larger in chronic liver disease9-13 than in normal liver13-16, which enables positive selection to shape the genomic landscape9-13. Here we analysed somatic mutations from 1,590 genomes across 34 liver samples, including healthy controls, alcohol-related liver disease and non-alcoholic fatty liver disease. Seven of the 29 patients with liver disease had mutations in FOXO1, the major transcription factor in insulin signalling. These mutations affected a single hotspot within the gene, impairing the insulin-mediated nuclear export of FOXO1. Notably, six of the seven patients with FOXO1S22W hotspot mutations showed convergent evolution, with variants acquired independently by up to nine distinct hepatocyte clones per patient. CIDEB, which regulates lipid droplet metabolism in hepatocytes17-19, and GPAM, which produces storage triacylglycerol from free fatty acids20,21, also had a significant excess of mutations. We again observed frequent convergent evolution: up to fourteen independent clones per patient with CIDEB mutations and up to seven clones per patient with GPAM mutations. Mutations in metabolism genes were distributed across multiple anatomical segments of the liver, increased clone size and were seen in both alcohol-related liver disease and non-alcoholic fatty liver disease, but rarely in hepatocellular carcinoma. Master regulators of metabolic pathways are a frequent target of convergent somatic mutation in alcohol-related and non-alcoholic fatty liver disease.

The mutational landscape of human somatic and germline cells.

Over the course of an individual's lifetime, normal human cells accumulate mutations1. Here we compare the mutational landscape in 29 cell types from the soma and germline using multiple samples from the same individuals. Two ubiquitous mutational signatures, SBS1 and SBS5/40, accounted for the majority of acquired mutations in most cell types, but their absolute and relative contributions varied substantially. SBS18, which potentially reflects oxidative damage2, and several additional signatures attributed to exogenous and endogenous exposures contributed mutations to subsets of cell types. The rate of mutation was lowest in spermatogonia, the stem cells from which sperm are generated and from which most genetic variation in the human population is thought to originate. This was due to low rates of ubiquitous mutational processes and may be partially attributable to a low rate of cell division in basal spermatogonia. These results highlight similarities and differences in the maintenance of the germline and soma.